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ABSTRACT 

In  the  heterogeneous  flow  of  finely  divided  particles  and  gas 
in  a  rocket  nozzle,  the  assumptions  that  the  flow  enters  the  nozzle 
with  no  lags,  that  the  radiative  transport  of  energy  between  particles 
is  negligible,  and  that  there  are  no  temperature  gradients  within  the 
particles,  are  commonly  made.     In  this  paper,  each  of  these  as- 
sumptions is  investigated  in  detail,  and  they  are  shown  to  be  reason- 
able for  typical  rocket  nozzles. 


I.     INTRODUCTION 

In  recent  years,  it  has  become  common  to  add  metals  to  the 
propellant  in  rocket  engines  to  increase  the  chamber  temperature. 
Unfortunately,  most  of  these  additives  form  some  sort  of  a  solid 
particle  in  the  exhaust,and  consequently  can  cause  a  depreciation  in 
the  thrust  and  specific  Impulse  from  the  values  based  on  gaseous 
products.     That  is,  instead  of  an  expected  20  per  cent  increase  in  im- 
pulse,  say,  an  increase  of  12  per  cent  is  noted. 

A  very  common  mixture  is  a  propellant  of  ammonium  per- 
chlorate  oxidizer  and  a  polysulfide  fuel  with  aluminum  added.     For 
the  numerical  work  in  this  thesis,   this  propellant  is  assumed  with 
pure  aluminum  added  to  it. 

The  amount  of  metal  that  can  be  added  is  limited,   since  there 
is  a  requirement  for  the  fuel  to  act  as  a  binder  for  the  oxidizer  and 
metal.     Usually,  the  designer  will  add  as  much  metal  as  he  can;  thus 
a  method  is  needed  whereby  he  can  correlate  test  engine  performance 
with  the  predicted  performance. 

Such  a  method  is  presented  by  Rannie  in  reference  1.     In  ref- 
erence 1,  Rannie  used  a  one -dimensional  hydraulic  approximation, 
which  is  known  to  be  valid  for  homogeneous  flow,  for  the  heterogeneous 
flow  of  gas  and  solid  particles.    In  addition,  he  assumed  (i)  that  the 
flow  entered  the  convergent  section  of  the  nozzle  with  no  lag  between 
the  gas  and  solid  particles;  (ii)  that  the  radiative  transport  of  energy 
between  particles  was  negligible;  and  (iii)  that  the  thermal  conductivity 
of  the  solid  particles  was  sufficiently  high  compared  to  that  of  the  gas 
that  there  were  no  temperature  gradients  within  the  particles. 
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As sumption  (i)  is  Investigated  in  Part  II  by  solving  the  equa- 
tions of  motion  in  a  cylindrical  combustion  chamber  and  determining 
the  velocity  lags  at  the  entrance  to  the  nozzle.     Axially  symmetric 
flow  is  assumed  in  the  chamber,  and  the  mixture  is  approximated  by 
a  heterogeneous  mixture  of  gas  and  solid  spherical  particles.     A 
constant  source  distribution  of  gas  and  particles  is  assumed  along 
the  cylindrical  wall,  and  viscous  effects  in  the  main  flow  are  ne- 
glected.    Turbulence  and  temperature  gradients  in  the  chamber  are 
also  neglected.     It  is  shown  that  the  particle  lag  at  the  entrance  to  the 
nozzle  is  indeed  very  small  compared  to  the  gas  velocity,   and  thus 
the  assumption  that  the  lag  is  zero  is  a  valid  one. 

The  magnitude  of  the  effect  of  radiative  transport  of  energy 
between  particles  is  estimated  by  adding  a  radiation  term  to  the  equa- 
tion for  the  heat  balance  of  a  particle  and  investigating  the  size  of 
this  additional  term.     The  ratio  of  the  absorption  cross  section  to  the 
geometric  cross  section  is  a  necessary  parameter  in  this  term,   and 
a  value  is  not  available  for  aluminum  oxide  particles  of  these  small 
sizes  and  high  temperatures.     However,  for  any  reasonable  value  of 
this  parameter,  the  radiative  transport  of  energy  is  shown  to  be  neg- 
ligible compared  to  the  heat  transfer  by  convection  and  conduction. 

The  temperature  distribution  within  a  spherical  particle  is 
found  for  unsteady  heat  flux  over  the  surface  of  the  particle  in  terms 
of  the  Laplace  transform.     By  expansion  of  the  transform,   the  solu- 
tion for  uniform  temperature  within  the  particle  is  recovered  and  the 
magnitude  of  the  correction  term  is  determined.     The  correction  term 
is  shown  to  be  negligible  for  practical  cases. 
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II.     TWO- DIMENSIONAL  ANALYSIS  OF  THE  FLOW 
IN  A  COMBUSTION  CHAMBER 

Following  the  notation  of  reference  1,  let    oC  be  the  mass 
fraction  of  particles  in  the  heterogeneous  mixture,     fa      the  gas  den- 
sity,  and     fa      the  density  of  the  solid  particles.     Then  the  density  of 

the  mixture  is  given  by 

_L         I  -<*r  <* 

r   =  T~      IT  <1' 

Let   A     be  the  cross-sectional  area  of  the  cylindrical  combus- 
tion chamber  and   m    the  constant  mass  flow  rate  of  the  mixture.     Ap- 
plying the  condition  that  the  net  flow  rate  into  any  volume  must  be 
zero,  the  continuity  equationsfor  the  solids  and  gases  separately  are: 

_cL_iV, api  i  1       i     I     h. 


oi 


[fl-*)nV|  ^~h\[[^)?^r\-  °  <2) 


l*f«s]     +  ^[*^ 


=  o 

(3) 


where     Ua     and     Us     are  the  axial  velocity  components  of  the  gas  and 
solids  respectively,   and    Vo      and     \ls      are  the  radial  velocity  com- 
ponents of  the  gas  and  solids. 

By  considering  the  momentum  transport  through  a  small  ele- 
ment and  neglecting  friction  forces,   the  momentum  equation  in  the 
axial  direction  can  be  written 


d 


x 


Zrl  J  OK 
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and  the  radial  momentum  equation  is 

tfr[(.-)f^v]  --^  +  f     (5) 

Approximating  the  particles  by  spheres  of  radius  cl    ,   and 
assuming  that  a  modified  Stokes  drag  law  holds  with  a  correction  fac- 
tor  fd        to  allow  for  deviations  from  Stokes'  flow,  the  force  balance 
in  the  axial  direction  on  a  particle  is 

^^■Tt^^h  ^M-^i  -  ^k        ,6, 

The  steady  state  drag  law  is  assumed  valid.   jU      ia  the  viscosity  of 
the  gas  and  is  constant  since  the  temperature  is  assumed  constant, 
and  the  force  due  to  the  pressure  gradient  is  added.     The  coefficient 
fj       is  a  function  of  the  Reynolds  number  and  the  Mach  number  of  the 
gas  relative  to  the  particle,   and  is  equal  to  1.  0  in  the  Stokes  flow  re- 
gime. 

A  similar  equation  for  the  radial  force  on  a  particle  is 

Equations  6  and  7  are  expressed  more  conveniently  in  the  fol- 
lowing form: 

^  Us~  <?  /*  r^^v  p,  *m  (8a) 
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and 


vf-vs 


-  iL  e^J-A  \  u.,iA  +  vf  ii  4-  J-  ids.) 


(8b) 


Letting    U5  =  U^  +  (U.s  -  U^)     and    l/s  —  \L  4- ( V5 - V-        in  equation 
3  and  subtracting  it  from  equation  2 


>S]+^[fvr]-lr[^^i:i;k(ur^f 


(9) 


Introducing  the  following  non-dimensional  variables  Y^-tt-    » 


"v^ 


r= 


_  j/. 


V0 


.  *  = 


_    -X_ 


b 


ii  =*   -£-  ,   where     V0       is  the  radial 


ve- 


locity component  of  the  gas  at   -n-  |      just  outside  of  the  boundary  lay« 
er  and    b    is  the  radius  of  the  burning  surface.     To  simplify  the 
equations,   also  let 

Introducing  these  in  the  above  equations,   we  obtain  for  the  combined 
continuity  equation 


3?Khf^ 


nt^f^-vM-^' 


-  o 


(10) 


and  for  the  momentum  equations 


(a) 
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and  for  the  force  balance  equations, 


-yjpi^^  ,12) 


41,  -  4c  -     e  J  4i 


H 


(13) 


(14) 


With  the  help  of  equations  2  and  3,   equations  1 1  and  1 2  can  be 
further  simplified  by  rearranging  so  the  dependent  variables  are    <£ 
\  i   T5  -  ta      »     and    f^-fa         .     Thus  equations  11  and  12  become 


■^ 


*F 


+  *,  *JdbL  +  (ft-^^-r,)^  =~it  <15) 


^^-  ^S^fa  ^^-W^p^1 


+ 


^^3l^-^  +  ^1)T^ 


^Obsfi  =  -±  4^- 


t  ^ 


? 


i 


\'o      ^ 


(16) 


First  Approximation,    Zero  Velocity  Lags. 

When  the  velocity  of  the  particles  approaches  that  of  the  gas. 


(i.e.,   as  (X.  — ^  o      ),    e.  — >  o       and    <p^ — =*- 4>o        an<*      % 


* 


* 
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Then 

Cps=  $  -  £0  and       rs=  ^  -  H,  (17) 

where  the  subscript  zero  Indicates  no  lag. 
Then  equation  10  reduces  to 

2>f  ^  ^  (18) 

and  equations  15  and  16  reduce  to 

^if  +  ^J  --fclf  (19) 


& 


(20) 


It  is  a  simple  matter  to  verify  that  the  solution  to  these  equa- 
tions,   satisfying  the  boundary  conditions    <^>    =  o  at  i  =  ° 

and     f0=-  —  /  at     t?=      I  ,   is 

Second  Approximation,    Small  Velocity  Lags. 

We  shall  assume  for  the  solution  that  the  variable  may  be 
expanded  in  powers  of    6      ;  thus 


<8- 


4>Q  =     +a  +  « 


J 


4>*  +-  <^  £v,   + 


=      oL0    \    6     oC,     + 


^ 


(22) 


; 


The  criterion  for  these  expansions  to  be  valid  is  that 


I 


(23) 


and  similar  inequalities  for  the  other  variables.     Here    6-    is  merely 
an  indicator  of  a  small  quantity.     Substituting  these  expansions  into 
equations  13  and  14,   and  equating  coefficients  of    £     , 


and 


^,-^. 


The  density  of  the  mixture,     f0       ,   is  very  much  less  than  the  density 
(%       of  the  solid  material,   so  those  terms  can  be  neglected  that  con- 
tain     ^°/ f<       *  therefore 


t 


h 


<fc.   =   4> 


* 


*   ^ 


"? 


(24) 


and  similarly 


-9- 


u>     _  V     _   4>   iJ£ 
T?.      7s,  -   ^    w 


^77 


(25) 


Now  defining  the  first  two  terms  of  equation  10  as    p    ,  then 


F  = 


Lf^]  +  H[f^] 


Equation  2  can  be  written  as 


f  -  it-T  * 


r^J  -f  A^rr^]  ~ 


Then  differentiating  this  by  parts, 

F  - 


*f*M-f#fw-f^-rr»^ 


=  o 


or 


BoL 


2) 


=   o 


F    - 


Also  comparing  equations  26  and  10  we  see  that 


^[<r(^]+^[ocr(v»s 


vi 


(26) 


(27) 


(28) 


From  the  first  approximation,   we  may  evaluate  equations  24 
and  25  as 


Ij,-  ^  =  >? 


(29) 


-10- 
Substituting  these  into  equation  28,   the  solution  of    F     is 

F=    l°cL0?0  (30) 

where  use  was  made  of  equation  22  and  only  carried  to  order  epsilon. 

Putting  equation  30  back  into  equation  27  and  incorporating 
equation  22,   and  rearranging. 

Particular  solutions  for  equation  31  are 

(32b) 

The  boundary  condition  for      <*-  ,      is  chosen  so  the. mass  flow 
rates  of  gas  and  solids  are  the  same  as  with  no  lag.     Since 


4 

A*.  =        <*- 


&*[,*.(**£+«»      '  <33) 


the  conditions  that     ^  3//4        and    ^f^A         are  tJie  same  as  with  no 

lag'are      rv=o 

sl±  +  a.  -t-  n    =  ©  (34) 

From  equation  1,  upon  substitution  of  equation  22,   and  dropping 

2- 

term3  of  order     €      , 
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P=    Pe>oC,  (35) 

Substituting  equations  35  and  29  into  equation  34,   we  find  that 

cLf^     -oLo[l-<to)  aX     y=   J  (36) 

satisfies  the  condition  that  the  last  term  of  equation  34  vanish. 
Clearly  the  solution  32a  is  not  applicable,  and  applying  condition  36 
to  equation  32b, 

DC  ,   =  -  L  oLo  (l-oLo)[Li  y>  f  -*-]  (37) 

The  boundary  layer  and  friction  forces  along  the  wall  are  ne- 
glected and  the  mass  fraction  of  particles  to  the  mixture  is  assumed 
to  be  constant  along  the  length  of  the  chamber  since  the  propellant  is 
presumably  uniform  in  mixture. 

This  solution  does  become  infinite  at  the  centerline,   so  the 
solution  is  valid  for  the  whole  chamber  except  on  a  line  down  the  cen- 
ter of  the  chamber.     This  singularity  does  not  cause  any  difficulty 
because  the  mass  flow  rate  in  the  axial  direction  is  finite,   as  will  be 
shown  below. 

Now  substitute  equation  2  2  into  equation  15,  and  put  -P^fo  +  Gf,* 
then  equate  the  coefficients  of  &     ,  dropping  known  zero  derivatives 
(i.e..    ^°/hti      and  ^/<3|         );  the  result  is 
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and  upon  substituting  equations  21,    29,   and  35, 

£fe  ?  +  ^  ^  >>  "  1^  - "  <*  -  ^ 

Combining  the  second  and  third  terms  and  putting  in  the  value  (equa- 
tion 37)  of    <*-/      ,  we  obtain 

-^^(^^;+4?fF^j-Y^-'^?=^4t   (3S) 

Similarly,   equation  16  becomes 

Let  us  assume  a  solution  of  equation  39  to  be 

Substituting  equation  35  into  the  continuity  equation  for  the  gas  (equa- 
tion 2)  leads  to 

Upon  substitution  of      ' ^/        into  equation  41  and  integration, 

f^    =     Ufl  f  X^    i-3°Lc  £  (42) 

where  the  constant  of  integration  disappears  since     <£<?,      must  be 
zero  at     f  —   ° 

From  equation  29  we  find 

and 


-13- 
To  verify  the  solutions,   substitute  the  values  for     (V     ,    T^/    , 
^       ,   and     Ts,        (equations  40,   42,   and  43)  into  the  continuity  equa- 
tions (equations  2  and  3).     The  mass  flow  rate  of  gas  for  a  cylindrical 
chamber  of  radius    b     ,  with  no  lags,  is 

fr  a    =       ft  if  ( I  -<o)  fo  Vo  <fe  (44) 

and  with  velocity  lags,  is 


=  x7U20-^Po^tJ(i^4!x)^7 


Substituting  from  equation  42  and  integrating,   equation  45  becomes 
m^    =    ^f  •  (46) 

The  mass  flow  rate  of  the  particles  without  a  velocity  lag  is 

and  with  a  velocity  lag, 


A^-  xtrb^ftv.*^'[l  +  fi(^i-§  +  ^]ilrfT 


(48) 


Evaluating  the  integral 


C5 


^   =>    m3  (49) 

showing  that  the  mass  flow  rate  of  the  particles  does  not  change  when 
velocity  lags  are  introduced.     Hence,  the  assumed  solution  and  con- 
stants are  consistent  with  mass  conservation. 

To  obtain  some  numerical  results,   a  specific  chamber  must 
be  chosen.     The  following  values  are  consistent  with  the  numerical 
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example  given  in  reference  1. 

b    =0.  3U      ft  CL^ZS^iO        cm. 

The  chamber  length  is  chosen  as  twice  the  diameter,   and  Stokes 
flow  will  be  assumed  (  J^-  /.o)  .     Then 

V6=     -12,7    fee  4>0  =    7^.  ?%ec 

6    =     3.H3  *  'e>  Yo  =  -I  9.7    /j€c. 

The  magnitude  of  the  velocity  lags  in  relation  to  the  no-lag 
flow  velocity  is  obtained  from  equations  29  and  21 

-^—^    =.i_e  (50) 

which  verifies  that  the  lags  are  indeed  small. 

The  uniformity  of  the  flow  across  the  exit  plane  was  assumed 
in  the  analysis  in  reference  1,  and  Table  1  is  presented  as  a  verifica- 
tion of  this  assumption.      Values  of    <=       fj^p^         ,    e       /  4  o         »   ^d 
^^'/cCc         were  calculated  as  functions  of  y     ,  the  non-dime  nsional 
radius. 

The  Reynolds  number  of  the  gas  relative  to  the  particles  was 
computed  at  the  exit  plane  using  the  particle  diameter  as  the  charac- 
teristic length  and  found  to  be  6.  5  X  10*   .     The  Mach  number  was 

-4 
also  found  to  be  1.  36x10     .     Combining  these  two  values,   a  value 

for    M/Re    of  2.  09  X  1Q~     is  found,   which  indicates  that  the  particles 

are  in  the  Stokes  flow  regime,   and  thus  the  correction  factor   fy        is 
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TABLE  1. 


? 

to 

•* 

6 

at, 

.  1 

-7.35  xio"4 

-3.69  x 

IO"4 

2.63 

X  IO"2 

.  2 

-6.  10  X 

io-4 

-2.  43  x 

io-4 

1.73 

X  IO"2 

.3 

-4.  75  X 

io-4 

-1.  55  x 

IO"4 

1.  28 

X  IO"2 

.4 

-4.  45  X 

io-4 

-9.  15  X 

IO'5 

9.  25 

x  IO"3 

.5 

-4.  07  x 

IO'4 

-4.21  x 

IO"5 

6.49 

X  IO"3 

.6 

-3.69  x 

IO"4 

-2.29  X 

io"6 

4.24 

x  IO"3 

.7 

-3.35  x 

IO'4 

+3.  14  x 

io-5 

2.  35 

X  IO"3 

.8 

-3.05  X 

io"4 

6.  15  x 

IO"5 

6.95 

-4 
x  10  * 

.9 

-2.31  X 

IO"4 

8.41  x 

io-5 

-7.  58 

-4 
x  10 

1.0 

-2.55  x 

IO"4 

1.  10  x 

io"4 

-2.06 

-4 

x  10  * 

approximately  equal  to  one,  as  taken  above.  Reference  1  gives  a  dis- 
cussion on  the  calculation  of  the  parameter  r ^  .  Using  equation  66 
and  equation  67  of  reference  1, 

L   ■=.     |.e>   -f   *SS  -**-  =     l.ooU  (51) 

Thus,  from  the  above  analysis,   the  assumption  made  in  refer- 
ence 1,   i.  e.  ,   that  at  the  entrance  to  the  nozzle  the  flow  was  uniform 
and  without  lags,   is  a  reasonable  assumption. 
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III.     RADIATION  BETWEEN  PARTICLES 

For  a  one-dimen3ional  steady  flow  of  a  homogeneous  fluid 
without  heat  addition,   the  energy  equation,   in  terms  of  internal  energy 
per  unit  mass,  is 

PfiU(e   t^    +     Afll    =   CentU,  (52) 

and  for  the  heterogeneous  mixture  of  gas  and  particles  is 

(i^m^^y^m^r^yAf  (ufo^f  ru,^)-  ««r.     |53) 

a 

where  oC  is  the  mass  fraction  of  particles  (subscript  s  ).  This 
equation  is  general  and  independent  of  force  interactions  and  heat 
transfer  processes  between  the  particles  and  the  gas. 

If  radiative  transfer  of  heat  between  particles  is  allowed,  an 
additional  term  must  be  added  on  the  left  hand  side  of  the  equation, 
since  the  equation  does  not  take  care  of  interactions  and  heat  transfer 
processes  between  particles.     The  energy  transfer  between  particles 
due  to  radiation  is 


whe  re 


k    =     Lk  (t~  o    j3 


3 


s  (55) 


may  be  regarded  as  an  effective  conductivity  due  to  radiant  energy 
transfer.     The  coefficient      0^      is  the  Stefan- Boltzmann  constant, 
and    Xn     is  the  Rosseland  mean  free  path.     Then    %r~~*'-(j      A 
must  be  added  to  the  left  hand  side  of  equation  53  to  account  for  radi- 
ative transfer. 


-17- 

Equations  54  and  55  are  applicable  only  when  -Zfz.      *8  very 
small  compared  to  the  length  scale  appropriate  to  the  problem,   for 
instance  the  nozzle  throat  radius,   at  all  important  wavelengths  of 
radiation. 

The  variable  mass  fraction    c<_    in  equation  53  can  be  elimi- 
nated with 

*  a,  =  (i  -c)  u9  f^o  (56) 

obtained  from  equations  2  and  3.     In  cases  of  interest,   the  volume  oc- 
cupied by  the  particles  is  negligible  when  compared  to  the  volume  of 
the  gas, and  the  assumption 

\~     0-4P  ,57) 

3 
is  valid  for  a  particle  material  density  of  210  lb/ft    as  we  have  chosen. 

3 
Since  the  gas  density  will  seldom  be  greater  than  0.  4  lb /ft  ,   this  ap- 
proximation is  valid  over  a  wide  range,   but  should  be  checked  for  any 
particular  problem. 

Substituting  equations  56  and  57  into  equation  53,   the  energy 
equation  becomes 


f1 

where    nc      is  the  enthalpy  in  the  chamber. 

For  a  perfect  gas  with  constant  specific  heats 


(58) 


y-f~  C^tRT^  =CpT§  (59) 

d 
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and  we  now  introduce  a  specific  heat,    6  Pc?     ,  for  the  mixture  with  a 
particle  mass  fraction     o^0    , 

■      ^?c  sa'    ^'-^)^p    +  *■*  ^  (60) 

and  modify  the  gas  constant  to 

R0=  ('-*•)  R  (6i) 

The  right  hand  side  of  equation  58  is     Lp0  1^       where    T^      is 
the  stagnation  temperature  of  the  mixture.     Then  equation  58  may  be 
rewritten  as 


tf-O  Vfc+okCTc  t  ^/j-^Vl^U^-ii    //_<,U2  =  ^po"k  (62) 

where  equations  56  and  59  were  used. 

Now  writing    Ts      as    U  t  1$  — T a        and    ^b        as    ^9  +  ^5-^3, 
and  combining  the  terms  of  equation  62, 


^  <"£  +*.4T«-T,>t  f(  if-itf)  -±rO^)g  =  ^ctc        (63) 

Neglecting  radiation,   the  heat  balance  for  a  sphere  of  radius 
<*-    and  temperature      1$      is 

f-Wft  C  U,  4g  =  -^  (Ts_T^  4  ^  (64) 

where     ka       is  the  conductivity  of  the  gas  surrounding  the  sphere  and 
j-\^     is  a  correction  factor  that  takes  into  account  the  influence  of  con- 
vection.    For  the  Stokes  regime  of  flow,  -f-,    =  \  .£>      ,   and  it  decreases 
at  higher  Reynolds  numbers. 
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Rearranging  equation  64, 


V^F^"^  (65, 


To  express  equation  65  in  dimensionless  form,   the  following 
notation  from  reference  1  is  introduced 


r  =     Tc   r 

€    _      i     a1^  f&7c~ 

1  rr/^ 


(66) 


giving 


where 

z9 =  £:  ^  K  <68> 

has  a  numerical  value  not  much  different  from  unity.     If  the  particle 
radius  is  extremely  small,    £     approaches  zero  and  since  the  other 
factors  on  the  right  hand  side  of  equation  67  are  of  order  unity, 
%~t%    —  °  a**d  the  temperature  lag  is  zero.     The  velocity  lag  is 

also  proportional  to   e      so  it  becomes  zero  as  well.     As  &-    increases, 
e      becomes  larger,  and  the  principle  of  the  perturbation  procedure  in 
reference  1  was  to  expand     '<?    ,      /j     ,    t  <l     ,   and  Ts    in  power 
series  in    £     . 
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V/ith  radiative  transfer,   there  is  an  additional  term  in  the  lie  at 
balance  for  the  particle.     Suppose  a  particular  particle  is  located  at 
X=  \0        .     The  radiative  transport  of  energy  in  the   X  -  direction  per 
unit  area  normal  to    X     in  the  space  durrounding  the  particle  and  close 
to  the  particle  is 


Hence  the  net  radiative  transport  of  energy  to  the  sphere  is 


1V 


»«*-«J««"'[-(kg^(ltM(<L,ta.)J 


-  -i^liC-Bi 


and  the  heat  balance  for  the  particle  becomes 


lW G  £C4d£  =  -  k(Ts-T*  WX  -^tr*  ^  (\<J& 


3J         t      3"      d- 


(69) 


which  replaces  equation  64.     Rearranging  and  introducing  the  dimen- 
sionless  variables  from  equation  66 


The  radiative  mean  free  path    i-  £  is  given  by 

4e,  L_ 


(71) 
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where    n.     is  the  number  of  particles  per  unit  volume  and    GL   is  the 
ratio  of  absorption  cross  section  to  geometrical  cross  section.     -Trom 
the  identity 

n-l-ft-^f,   =   ^f3  ,72, 

and  the  equation  of  state  for  the  gas 

nfT*=3-SL.   _£- .      =    3.     ■& «     T^  (73) 

Substituting  into  equation  71  above 

^=£«SP^-i-    ^f^  (74) 

and  hence      kK     from  equation  55  becomes 


Substituting  for     ^        into  equation  70 


Vn--4>^4Jl^^ft%      fitfe'f^^fl    (76, 

where    Q    is  assumed  to  have  small  variation  with     f     and  can  be 
taken  out  of  the  differential.     Since    ^f/fc    -    VX~(  »   equation  76 

can  be  written  as 


1  "*  "  *  ^  £  ^  $  ^  p^  ^(¥  ~^if)j 


(77) 


where  r  _    JlSL    *-_<_ 


-22- 

Slnce  the  factors  multiplying     j-       in  equation  77  can  be  con- 
sidered to  be  of  order  unity,   ae  is  the  first  term  in  the  square  brackets, 
a  value  of     4>      o£  order  unity  implies  that  radiative  transport  is  as  ef- 
fective as  the  heat  transfer  by  convection  and  conduction.     In  the  exam- 
pie  of  reference  1,  -f>^  3*1°/  Q.     ;  hence,    yv     is  extremely  small 
even  if    GL  is  quite  small.     Hence  the  influence  of  radiative  transport 
on  particle  lags  is  negligible. 

The  zero  lag  solution  is  also  affected  by  radiative  transport  of 
energy.     The  energy  equation  63  becomes 

Upon  substitution  of  equation  55,  we  find  that  the  condition  for  radiation 
to  have  negligible  effect  on  the  zero  lag  solution  is 

^'^       Jfi      ^^     j  (78) 


-3 


since  the  non-dimensional  terms  are  all  approximately  of  order  one. 

The  concept  of  an  effective  conductivity  has  meaning  only  if 
^V/-r    **    I         •     For  chemical  rockets    ^VtfjTeli)  =     /<? 
Hence  the  radiative  term  gives  negligible  contribution.     The  ratio 

In  the  example      <^/r  -S!   3  *   /  o~ "      ,    '?  ^    ■   £*  5T  fi  1°       ,   and 


/-*    =z    US'     •     Sven  if     &•  is  small,     -%/^    ^-C  / 


OV 
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IV.     EFFECT  OF  PARTICLE  CONDUCTIVITY 

The  effect  of  finite  particle  conductivity  is  investigated  by 
solving  the  heat  conduction  equation  for  the  particle  with  an  appropri- 
ate boundary  condition,   then  expanding  that  solution  to  investigate  the 
size  of  the  correction  term. 

The  heat  conduction  equation  for  a  spherical, uniform, solid 
particle  with  spherically  symmetric  temperature  distribution  is 


AT  =  a-  vY-  cr-L  4-  (  r^ 


(79) 


where   CT    is  the  diffusivity  of  the  particle  material. 
The  initial  condition  is 

T=  Ta  at     t  =  o  (80) 

and  the  boundary  condition  at  the  particle  surface  is 


k'(£Uu.-=fc»(#) 


by-jr-o^-o  jv^ir-ato  (81) 

where       ic.     is  the  chamber  temperature  and  ~£    is  the  time  measured 
from  the  instant  the  particles  enter  the  convergent  portion  of  the  noz- 
zle.    The  conductivities  of  the  solid  and  gas  are       j^.s     and      lc*     re- 
pectively. 

For  steady  6tate,   with  spherical  symmetry,   the  equation  for 
heat  conduction  in  the  gas  is 

^v-l        *y  J         °  (82) 
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which  has  the  solution 

T=   -A-+-  B 

r 

At    1^-    °o       ,    T==  Ta  ,   and  thus    B  =  Tq       and  at    IT"— 

T=    ~T*_    which  gives  A  =  cl  .( "^-~" "^ )  •     Thus 

(ft)       --i(Trt) 


CL- 


OT 


k»^La*.-      %^T-) 


(83) 


For  quasi- steady  heat  flow  in  particle  and  gas,    substitute 
from  equation  83   into  81     to  obtain 


This  will  be  used  for  the  approximate  boundary  condition  for  the  heat 
conduction  equation  of  the  particle.     The  solution  to  these  equations 
(79  t  80     f  84     j  -wiH  be  obtained  by  making  use  of  the  Laplace  trans- 
form. 

Define 

f  \ef  pCOJt  =Z[HcJ\=   TZf)  (85) 

and 

rT~    —     ®  (36) 
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Then  equation  79   becomes 


&JI-    fr  &   +  J|T   Tc—  o  (87) 

3T" 
Since  there  is  no  source  at    T"  —O       ,  >— .  =  o      at    V  =  °       and 

therefore   ~T  -  °  where     T  ~   c  ~VZ  r  =  ~  *^>  +"  G4~ 

which  means  that 

r  &  -  3  — >  <=>  as     y =>  o  (88) 

The  solution  of  equation  87  is 


o  =  4 e  t    Be  +  T c  r  (89) 


Substituting  this  into  equation  88 


he      -t  Be  +  ToW  1VU' 


and  to  satisfy  condition    83    ,     D  —      A     t  thus 

T«  ■£-"  =  A  a.  sink  a*-  +  Ta 

where  a  =   j^O 

Substituting  this  into  the  boundary  condition  (equation    84  ), 
we  obtain  for  the  left  hand  side 


(91) 


y-  -c3^ 
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ajtid  for  the  right  hand  side 


rf-^j-T*.)  -  *|t, 


=      U  r  T.  -  ^Alinlvyo-  __-p  "I 


<a- 


since    (     =    To*,         at    ^-a.      .     Combining  these  equations,   equa- 
tion 84    becomes 


44-H^- ~si^]l  -  -^f  Tg  -  2=^F  -Tc" 


(92) 


and  solving  for  A 


A  =    cu  k.(T9-Tc) 


34-^L 


*s  r  ^  d^s/,  aa-SiflL^  [[  -JS5- 1 


(93) 


Substituting  equation  93    into  91     gives 


t=«4(t^' 


y-k 


rink 


* 


^flLtojk«flL-5«n/i  ««-£<•  --=i\ 


+ 


(94) 


Let  us  now  define  an  average  sphere  temperature     15    as 


X 


J"    tlrSTAr 


and  hence 


^=mVt^ 


(95) 
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Then     ~TS      will  be  compared  with  the  uniform  sphere  temperature  in 
the  heat  balance  equation  for    lc;s      infinite  (reference  1). 


i^S^-t^WaJjL  ,96, 


Also  define: 


B=      *.-£* (T^TC 


<5     Qa-coskv^-  Sir\kaa.(l  -  ^^j 
so  that  equation  94  may  be  written  as 


T  =     8  £1^4*    +    Tc  (97) 


Then  substituting  and  evaluating  equation  95 


Ts     =     3   B,j    cosIi^a.    -   S-inK  fr*A    4-     7^,  (98) 


Because  large  values  of  time    "C    correspond  to  small  values 
of  <2-    ,   expand  the  egression  for       lj      in  a  power  series  in  <2_     to 
obtain  the  form  correct  for  small    a^  ,  i.  e. , 

ix-T^\^ui»r+kH*^Tik-M+ „fT    (99) 
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Retaining  only  terms  to  order     (<2  a.)      in  equation  99 


which  reduces  to 


rrT4  =  (Ts-r0|it^i+^a;_J_(?a/J 

Recalling  that  a^  r£r   =         ,  we  find  that  for 

equal  to  infinity,  equation  100  reduces  to 


(100) 


(101) 

and  talcing  the  inverse  Laplace  transformation  of  this  leads  to 

^l     1  J        3  ^  (10?) 

which  is  identical  with  the  assumption  used  in  the  heat  balance  of  a 
particle  of  infinite  thermal  conductivity,  i.e.,  equation  96. 

Reference  2  gives  values  for  the  thermal  conductivity  of 
aluminum  oxide  and  its  variation  with  temperature.     Measured  values 
are  given  for  temperatures  up  to  1600   C,   and  one  more  point  (1300   C) 
is  presented  from  extrapolated  data.     The  thermal  conductivity  de- 
creases with  an  increase  in  temperature  up  to  1400   C,   and  then  be- 
gins to  increase  again.     The  value  given  for  1400   C  (the  minimum 

value  of      ^5    for  this  range)  is  .0125 * .  ~ —    or 

cm     sec  C 
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36.25*™.  J2.     . 
fThr       °F 

To  find  a  reasonable  value  for  the  thermal  conductivity  of  the 

gas,   the  data  from  the  numerical  example  in  reference  1  was  used, 

that  is, 

C      =    0.500  BTU/lb  °F, 
P 

n      =    1.  5  x  10"6  lb.  sec/ft2, 

Pr    *    0.74    . 
To  be  most  conservative  in  evaluating       s/k.a       $   a  value  of 
yu_  was  used  that  would  give  a  high     k-«       .     Since  ft.     varies  con- 
siderably with  temperature  and  the  Prandtl  number  is  almost  constant 
over  a  large  temperature  range,   the  thermal  conductivity  varies  with 
the  temperature  in  the  same  manner  as  ^u.    .     Over  the  temperature 
range  of  interest,   a  power  law  of  the  following  form  is  satisfactory. 

^  «  Cs.  u-a  L5l\  (103) 


i  -  t^ 


The  maximum     Ic  <l     would  be  obtained  when    T^  =  lc_       »   thus 


Thus  the  value  of        %/ka        in  equation  *00  will  always  be 
greater  than  25.  7  ,   and  thus  one  may  neglect  1/5  in  comparison  to 

With  this  simplification  in  equation  *^°  and  substituting 


^K, 


for    ^  *"     j 
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3  i^-f    S4-    ^    r] 


Talcing  the  Inverse  of  the  Laplace  transform  yields 


(104) 


(105) 


Without  the  second  term  on  the  right,   this  equation  is  identical  with 
equation  101»   and  the  last  term  is  a  correction  term.     To  Investigate 
the  size  of  the  last  term,   rewrite  equation  102   a3 


and  differentiate  to  obtain 


&=at-i^-^  (107) 


_       3'^ 


where  oL  —  )/%(lcC'  .     Then  equation  105  becomes 


(108) 


The  condition  that  the  approximation  used  in  the  heat  balance 
equation  of  a  particle,   equation  96     t   is  valid  is 

_3    /Uf^-iB 


a-?  IwJ  Jd^         dr  '  ^<_    j  (109) 

It 
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For  a  first  approximation,   one  can  say  that         f'/af      and 

J  "7~ 

"      /&T        are  of  the  same  order,  while  still  unevaluated.     The  coef- 

•4 
ficient  of  the  numerator  is  1.  621  x  10      based  on  the  constants  evalu- 
ated earlier.     Since  ^  %Q,r  ~~       /<It    can  be  at  most  the  same  order 
as  Vdr         a-nd  will  probably  be  less,  we  can  say  that  condition  109 

is  satisfied. 

Since  this  condition  is  satisfied,  equation   96  is  a  valid  ap- 
proximation which  would  be  in  error  by  much  less  than  one  per  cent. 
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LIST  OF  SYMBOLS 

f\  cross  section  area 

d  specific  heat 

/^  Mach  number 

Pr  Prandtl  number 

^  ratio  of  absorption  cross- sectional  area/geometrical 

cross -sectional  area 

f\  gas  constant 

j^e  Reynolds  number 

I  temperature 

IL  axial  velocity  component 

V  radial  velocity  component 

(X-  radius  of  solid  particle 

P  radius  of  combustion  chamber 

£  internal  energy /unit  mass 

T  correction  factor 

la  enthalpy 

k  thermal  conductivity 

£n  Ross  eland  mean  free  path 

yh  mass  flow  rate 

n  number  of  particles /unit  volume 

"p  pressure;  variable  of  Laplace  transform 
rate  of  heat  transfer 
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% 

w 

r 

radial  distance;  burning  rate 

t 

time 

X 

axial  distance 

OC 

mass  fraction  of  particle 

t 

ratio  of  specific  heats 

e 

indicator  of  a  small  quantity 

r 

pressure  ratio 

\ 

non-dimensional  radial  distance 

& 

angle    ( radians ) 

H* 

viscosity 

f 

non-dimensional  axial  distance 

r 

density 

Cr" 

diffusivity 

« 

Stefan- Boltzmann  constant 

r 

temperature  ratio 

* 

non-dimensional  axial  velocity  component 

r 

non-dimensional  radial  velocity  component 

Subscripts 

J)  burning  surface 

q  chamber 

d  drag 

q  gas  phase 

h  convection 

£  with  lag 
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f 

constant  pressure 

t 

radiation 

s 

solid  particle 

t 

throat 

V 

constant  volume 

o 

no  lag 

1 

with  small  lag 

